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ABSTRACT 



Context. The wings of Balmer lines are often used as effective temperature diagnostics for late-type (F, G, K) stars under the assump- 
tion they form in local thermodynamic equilibrium (LTE). 

Aims. Our goal is to investigate the non-LTE excitation and ionisation of hydrogen and the formation of Balmer lines in late-type 
stellar atmospheres, to establish if the assumption of LTE is justified. Furthermore, we aim to determine which collision processes are 
important for the problem; in particular, the role of collision processes with hydrogen atoms is investigated. 

Methods. A model hydrogen atom for non-LTE calculations has been constructed accounting for various collision processes using 
the best available data from the literature. The processes included are inelastic collisions with electrons and hydrogen atoms, mutual 
neutralisation and Penning ionisation. Non-LTE calculations are performed using the MULTT code and the MACKKL semi-empirical 
solar model, and the relative importance of the collision processes is investigated. Similar calculations are performed for MARCS theo- 
retical models of other late-type stellar atmospheres. 

Results. Our calculations show electron collisions alone are not sufficient to establish LTE for the formation of Balmer line wings. 
Mutual neutralisation and Penning ionisation are found to be unimportant. The role of inelastic collisions with neutral hydrogen is 
unclear. The available data for these processes is of questionable quality, and different prescriptions for the rate coefficents give sig- 
nificantly different results for the Balmer line wings. 

Conclusions. Improved calculations or experimental data are needed for excitation and, particularly, ionisation of hydrogen atoms 
in low-lying states by hydrogen atom impact at near threshold energies. Until such data are available, the assumption of LTE for the 
formation of Balmer line wings in late-type stars is questionable. 

Key words, line: formation 



1. Introduction 

For more than 80 years it has been known that hydrogen is the 
dominant chemical element in the Sun, and most other stars, by 
number of atoms as well as by weight. For the Sun, and even 
more so for more metal-poor stars of solar type, hydrogen is the 
dominant continuous opacity source through H and H i bound- 
free and free-free absorption, as well as Rayleigh scattering. For 
the metal-poor stars of solar temperatures, hydrogen is also im- 
portant as an electron donor and thus governs the H opacity. In 
addition, the hydrogen lines in late-type stars are important as 
temperature diagnostics, and the Balmer discontinuity as a grav- 
ity diagnostic for F stars. 

When the development of electronic computers and suitable 
algorithms in the 1960's made it possible to reasonably self- 
consistently model excitation and ionisation together with radia- 
tive fields in stellar atmospheres, these advances were immedi- 
ately applied to the case of hydrogen. Most studies considered 
early-type stars, but the situation for late-type stars was also ex- 
plored (see Mihalas & Athay [1973] for a review of the early de- 
velopment). The newly developed Feautrier method was used by 
Cuny (|19671 l to model hydrogen in the solar atmosphere with a 
3-5 level model atom with continuum. The Ha line was found 
to be controlled by the photoionization processes from its lower 
and upper level, and to be sensitive to the structure of the photo- 



sphere. Strom (1967 ) studied the effects of departures from lo- 
cal thermodynamic equilibrium (LTE) for hydrogen on continua 
of G and K stars. The equations of statistical equilibrium were 
solved for the three lowest levels (assuming detailed balancing 
in the lines) and FT". Observable effects on fluxes and colours 
were found, but these vanished when more realistic values for 
the rate constant of the important associative detachment reac- 
tion, H" + H ^± H 2 + e, were used (Lambert & Pagel 1968; 
Pagel fl968l . 

Vernazza et al. (1981) calculated the statistical equilibrium 
of hydrogen in their semi-empirical solar models using a 12- 
level model atom with continuum. They included resonance and 
van der Waals broadening, but neglected Stark broadening, as 
did Cuny ( 1967). Departure coefficient^]/? significantly different 
from unity were found in the upper photosphere. When compar- 
ing their hydrogen line profiles to observations, they found rather 
severe discrepancies for Ly/3; for Ha, the results in the core were 
shown to be dependent on the number of levels used as well as 
on the particular model atmosphere. Vernazza et al. also included 
statistical equilibrium calculations for H in their solar models. 
They found departure coefficients that are significantly greater 
than unity for the continuum optical depths log Tsoonm < -5, and 
0.9 to 1.0 in the depth range -4 < logTsoo < _ 2, three times 
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1 In this paper we define the departure coefficient of level i as 6, = 
«,/«*, with n, the population and n* the LTE population. Some other 
authors use the definition b\ = (nj/n*)/(n K /n* K ), where k is the ionised 
state. 
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greater than those found by Lambert & Pagel ( 1968). The dif- 
ferences seem to be mainly caused by the fact that the tempera- 
ture minimum of the solar model atmosphere used by Vernazza 
et al. (U981I I (their model C for the average quiet Sun) is deeper 
and located further out than in the model used by Lambert & 
Pagel (|1968l l. At log rgoo > -2 the bu- coefficients found by 
Vernazza et al. (1 198 li t are unity, verifying the conclusion of 
Lambert & Pagel (1968 ) that LTE is a good approximation for 
the calculation of PL continua. 

Recent improvements in high-resolution stellar spectroscopy 
have made it possible to obtain observed stellar Balmer line 
profiles with reasonable accuracy, stimulating further study of 
these lines. The formation of the Balmer line wings occurs at 
relatively large optical depths compared to most of the spec- 
tral lines in late-type stars. This is due to the high excitation 
energy of the n — 2 level, and makes it possible to probe the up- 
per convection zone of the atmospheres. Fuhrmann et al. ( 1993 ) 
compared observed Balmer line profiles for late-type dwarfs 
with calculated profiles, based on the LTE approximation. They 
found discrepancies which were attributed to effects of inade- 
quate temperature structures, as obtained with a mixing-length 
description of the convective energy transport, and suggested 
revisions of the mixing-length parameters. In following papers 
(e.g. Fuhrmann et al. 119941 ), fitting of Balmer line wing profiles 
was used to obtain effective temperatures. A number of studies 
of Balmer line profiles in late-type stars in relation to convective 
parameters and temperature scales followed; for example, van't 
Veer- Menne ret et al. 019981 ), Gardiner et al. ( 119991 ), and Barklem 
et al. ( [20021 . 

All these recent studies, have assumed that the Balmer wings 
form in LTE. Vernazza et al. calculated departure coefficients of 
unity for all levels below logT5oo ~ -0.5; however, departures 
are seen in higher regions of the atmosphere. Similar results 
were found by Carlsson & Rutten (119921 ) in their study of so- 
lar Rydberg lines. Both of these studies included only collisions 
with electrons in the model atoms, and did not compare com- 
puted Balmer lines with those computed in LTE. As the Balmer 
line wings typically form in the region around -2 55 log T500 S= 
0.5, based on these results one would expect departures from 
LTE in Balmer line wings. The assumption that collisions with 
electrons are the dominant collision process is based on the fact 
that at the thermal velocities of interest, electron collisions are 
expected to be non-adiabatic and have substantial cross sections, 
while collisions with heavier particles, such as H, are expected to 
be nearly adiabatic and thus have small cross sections. Such ar- 
guments can be made from the Massey criterion (Massey 1949 
see also Anderson 119811 for an illustrative discussion), compar- 
ing the typical collision time with the inverse of the natural fre- 
quency of the transition (given by h/ AE, the period of light cor- 
responding to the transition). If the collision time is similar to 
this period then the collision is non-adiabatic and a transition 
is relatively likely. If the collision time is much larger than this 
period then the collision is nearly adiabatic and a transition is 
expected to be unlikely. Furthermore, the high thermal velocity 
of electrons leads to a higher collision rate. However, even if the 
inelastic cross sections for collisions with heavy particles such 
as hydrogen atoms are indeed small, the abundance of perturbers 
is extremely large. In the photosphere of solar-type stars neutral 
hydrogen atoms typically outnumber electrons by a factor of 1 4 , 
and can be orders of magnitude greater still in metal-poor stars. 

A recent study of the solar Balmer and Paschen lines by 
Przybilla & Butler (2004b ) did compare LTE and non-LTE line 
profiles and found them to be identical in the wings. While their 
model atom employed improved collision rates for electrons, it 



also included collisions due to hydrogen following the semi- 
empirical formulae of Drawin (119681 fl969) empirically scaled 
by a factor of two in order to reproduce the observed solar line 
profiles. 

Thus, it has never been established that Balmer line wings 
in late-type stellar spectra are formed in LTE. The basis of the 
assumption is that the line wings form deep in the atmosphere 
where collisional processes are assumed to be efficient and dom- 
inant over radiative processes. Even if this turns out to be the 
case, the atomic processes responsible have never been identi- 
fied. In this work, we make a non-LTE study of hydrogen and 
the formation of Balmer lines in late-type stellar atmospheres. 
In particular, we examine a number of collisional processes, es- 
pecially those involving perturbing hydrogen atoms, and their 
effects on the statistical equilibrium of hydrogen using available 
estimates of the collisional transition probabilities. Modern line 
opacity data are employed. In order to identify the important col- 
lision processes and examine the validity of the LTE approxi- 
mation in the calculation of Balmer lines for late-type stars, we 
have undertaken statistical equilibrium calculations in the semi- 
empirical solar model of Maltby et al. (1986) and MARCS theoret- 
ical model atmospheres of dwarfs and giants of different metal- 
licity (§[2]and[3](. Finally, we discuss the results and our conclu- 
sions (§|4]). 

2. Statistical-equilibrium calculations 

The equations of statistical equilibrium and radiative transfer 
are solved for a Hi model atom with 19 bound levels and con- 
tinuum in a ID plane-parallel model atmosphere, consistently 
with the hydrostatic equilibrium equation, using the MULTI code 
(Carlsson 1986 Carlsson et al. 119921 ). The atomic data for the 
model atom will be described in detail below. 

The calculations assume complete redistribution of the line 
radiation for hydrogen lines. Hubeny & Lites d 19951 ) studied the 
effects of partial redistribution on hydrogen in the Sun. Their 
method was also applied by Sim (2001) to other late-type stars, 
including Procyon (F5 IV- V), /? Gem (KO III) and a Tau (K5 
III). Partial redistribution is particularly important in strong res- 
onance lines, such as Lya. Although the direct effects on the line 
profile are most obvious, and lead to effects, for example, on 
the ionization balance, some more secondary effects also appear 
via the populations of excited states. Partial redistribution effects 
may also be significant via "cross-redistribution" (or resonant 
Raman scattering) between Ly/3 and Ha, sharing the same upper 
level. That is, a photon may scatter in one line and reappear in the 
other before the atom undergoes a redistributing collision. These 
effects were included in the calculations of Hubeny & Lites, and 
in those of Sim, but they do not include partial redistribution di- 
rectly in the subordinate lines like Ha, arguing that the effects of 
partial redistribution in these lines are negligible for the condi- 
tions in the solar atmosphere (see also Hubeny & Heinzel [T984| ). 

2.1. Model atom 

2.1.1. Energy levels and radiative transition data 

We employ a 19-level-plus-continuumHi model atom based on 
the model atom of Carlsson & Rutten (119921) . from which all 
energy levels and radiative transition probabilities are adopted. 
Numerical experiments suggest that results converge with in- 
creasing number of levels (e.g. Przybilla & Butler 2004a), and 
for solar- type stars 19 levels is found to be more than sufficient 
to achieve such convergence in line forming regions. 
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Detailed line profiles are required for the statistical equilib- 
rium and radiative transfer calculations. The hydrogen line pro- 
files are calculated including Stark broadening, self-broadening, 
fine structure, radiative broadening, and Doppler broadening 
(both thermal and turbulent). The line profiles of the Balmer 
lines of particular interest here (Ha, H/3 and Hy) are described 
as accurately as possible, including an accurate convolution of 
the profiles resulting from the different broadening mechanisms. 
The Stark broadening profiles from Stehle & Hutcheon (1999) 
are employed, but if the plasma parameters (T, N e ) exceed the 
range for which Stehle & Hutcheon have tabulated line profiles, 
the Stark broadening profiles from Vidal et al. (1973) are used 
instead. Self-broadening from Barklem et al. (2000) is used. The 
Gaussian (Doppler) and Lorentzian (self-broadening and radia- 
tive) components are convolved together employing the Voigt 
function, which is then convolved with the fine structure. The re- 
sulting profile is then convolved with the Stark broadened profile 
using the Fast Fourier Transform technique for computational 
speed. 

All other hydrogen line profiles are calculated in a more ap- 
proximate manner for fast computation. We use a slightly re- 
vised version of the code HLINOP which has been described 
by Barklem & Piskunov (2003), which is based on the original 
HLINOP by Peterson & KuruczQ. In this case, the Stark broad- 
ening is calculated using the theory of Griem (1960 and sub- 
sequent papers) with corrections based on Vidal et al. ( 11973I ). 
Self-broadening is included following the resonance broadening 
theory of Ali & Griem ( 1966). The profiles are convolved in an 
approximate manner. In the line core the Stark profile is approx- 
imated by a Lorentzian, and thus convolution is approximated 
via a Voigt profile. In the line wings, profiles are simply added 
together. We made test calculations comparing the results for Ha 
and H/? if computed in the detailed manner described above, with 
the results where the lines are computed with the more approx- 
imate method just described. The differences in the computed 
line wings for solar-type stars are very small, less than 0.08 % 
and 0.3 % in relative flux for Ha and H/3, respectively, in the 
solar case. The differences in the line cores are more substan- 
tial, particularly in the widths of the cores; at the sides of the 
cores the differences are of the order of 5 % and 2 % for Ha and 
H/3, respectively. Thus, a more careful convolution of profiles is 
needed for accurate calculation of the line cores. 

2.1 .2. Collisional transition data 

The collisional transition data has been significantly updated and 
expanded from the model atom of Carlsson & Rutten (1992 ), 
which considered only collisions with electrons. We added ad- 
ditional collisional processes in order to assess their importance. 
These processes and their modelling are detailed below. 

For some processes, rate coefficients could be taken directly 
from the literature, and in these cases the references are provided 
and may be consulted for the relevant data. However, in other 
cases we have derived rate coefficients from cross sections avail- 
able in the literature via numerical integration over a Maxwellian 
distribution of velocities. Where appropriate, we provide the rate 
coefficients via a fit to numerical results of the form 

log(crv) - a\+ 02 log T + fl3(log T) 2 + fl4(log T) 3 . (1) 

The fits are valid over the fitted temperature range of 1000 < 
T < 20000 K. 



2 see http://kurucz.harvard.edu/ 



3 

The collisional processes accounted for are now detailed. 
Rate coefficients for reverse processes are calculated from the 
principle of detailed balance. Throughout, n denotes the princi- 
pal quantum number of the level. 

Collisional excitation and ionisation by electrons 

H(n) + e ^ H(n') + e (2) 

H(n) + e ^ H + + 2e (3) 

Data from Przybilla & Butler (2004a), computed using the 
/^-matrix method in the close coupling approximation, were 
used for transitions between states with n < 7. For colli- 
sional ionisation, experimental results are available for ionisa- 
tion from the Is state from Shah et al. ( 119871 ), and from the 
2s state from Defrance et al. ( 119811 ), which are in good agree- 
ment with calculations, e.g., Bray & Stelbovics (1993) and 
Mukherjee et al. ([1989), respectively. We computed rate coef- 
ficients from the experimental cross sections; polynomial fits 
to the results are provided in Table Q] Calculations (see, e.g., 
Witthoeft et al. 2004) indicate that the cross sections for excita- 
tion from the 2p state are of the same form as for 2s, although 
approximately 20% larger. Thus, assuming the 2p rate coefficient 
to be a factor of 1 .2 larger, for the n — 2 state we adopt the 2s 
result multiplied by a factor of 1.15, which combines the 2s and 
2p results according to their statistical weights (2/8 +6/8 x 1.2). 

For excitation and ionisation involving higher states (n > 7 
for excitation, n > 2 for ionisation), formulae from Vriens & 
Smeets (1980) are used, as in Carlsson & Rutten ( 1 19921 ). These 
analytical formulae are a semi-empirical combination of experi- 
mental and theoretical data. Przybilla & Butler (2004a) made a 
comparison of available data for excitation and found the data 
of Percival & Richards ( 1978) to be the best alternative among 
several well-established approximations (not including Vriens & 
Smeets). Noting that for the case of hydrogen the expressions of 
Percival & Richards ( 1 19781 ) are equivalent to the expressions of 
Gee et al. ( 119761 ), Vriens & Smeets ( 1980) compare extensively 
with these results and find excellent agreement (see, e.g., Figs. 6- 
8 of Vriens & Smeets |T9801 ), and claim a wider range of validity. 
Figure Q] compares the ionisation rate coefficients from Vriens 
& Smeets ( 1980 ) with those derived from experimental data for 
the two lowest states discussed above. The agreement is very 
good at the temperatures of interest, giving confidence that this 
approximate formula probably gives reasonable results for more 
excited states. 

Collisional excitation and ionisation by H atoms 



H(«) + H(ls) ^± H(n') + H(ls) (4) 

^ H(ls) + H(n') (5) 

H(n) + H(ls) ^± H + +H(ls) + e (6) 

^± H(ls) + H + +e (7) 



The excitation process («' > «) includes two possible channels, 
the direct excitation process (HJ) and the excitation process in- 
volving excitation transfer ||5). Similarly, ionisation includes the 
direct process (O and the channel involving an electron trans- 
fer (Q. We note also the existence of the associative ionisation 
process H(n) + H(ls) ^ HJ + e, which is not included. 

Data for these processes are available from a number of 
sources. As we will see, comparison of the data from these 
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Fig. 1. Comparison of rate coefficients for ionisation of n = 1 
(upper panel) and n — 2 (lower panel) levels of hydrogen. The 
full lines are derived from the experimental results discussed in 
the text, and are used in the calculations. The dashed lines show 
the result from the approximation of Vriens & Smeets (11980I I. 



sources give a wide range of values for the rate coefficients, and 
it is not always clear which set of data are to be preferred in 
each instance. We will discuss each set of data in turn, and then 
compare the data for some example cases. 

Bates & Lewis ( I7955D : Data for the n — 3 — > 2 process have 
been calculated by Bates & Lewis ( 1955) using the Landau - 
Zener model. These calculations consider the ionic curve 
crossing mechanism, which is only important for this one 
excitation transition in hydrogen. This process may occur 
only between the n = 2,3,4 levels as only these states have 
avoided crossings with the ionic state potential; the n — 4 
crossing is at large internuclear distance and is thus ineffi- 
cient (Bates & Lewis [19551 ). A fit to the rate coefficient for 
the reverse process is given in TableQ] 

Drawin {1968 1969): General analytic expressions for the 
cross sections and rates for both excitation and ionisa- 
tion are given by Drawin ( 1 19681 119691 see also Drawin & 
Emard 119731 1. The expressions are derived from a semi- 
empirical modification of the classical Thomson formula for 
ionisation by electrons. Fleischmann & Dehmel ( \912\ have 
pointed out that the derivation has some apparent "inconsis- 
tency and arbitrariness". For example, despite the fact that 
Drawin states clearly that the atomic kinetic energy E a is in 
the centre-of-mass frame, this is inconsistent with the stated 
threshold of twice the ionisation energy which would in- 
stead imply that E a is in the laboratory frame. Furthermore, 
if we numerically integrate rate coefficients from Drawin's 
cross sections, we find that good agreement with his expres- 
sion for the rate coefficient can only be obtained if we as- 
sume E a is in the laboratory frame. Moreover, comparison 
with experiment for the case of ionisation (see, e.g., Fig. 3 
of Kunc & Soon 119911) . shows that while the formula may 
give reasonable results for collisions near the peak in the 
cross section at about E\. d \, ~ 10 4 eV, the results at lower 
energy (100-1000 eV) are not in good agreement. Note that 
for the temperatures of interest the near threshold collisions 
(£)ab - 2E cm = 27-40 eV) are of interest. 



Soon (1992): Similarly general analytical expressions for the 
excitation and ionisation cross sections and rates have been 
derived by Soon d 19921 1, that are an extension of the work of 
Kunc & Soon (1991) based on the classical impulse approx- 
imation (Gryzinski |19651 l. 

Unfortunately, we find problems with the presented formu- 
lation in this case also. The variable E is the impact energy, 
and appears to be defined as the kinetic energy in the labo- 
ratory frame, since it is stated that the relative kinetic energy 
in the centre-of-mass frame is given by E cm = E/2. If we 
take the example of the ionisation case, the cross sections 
have a threshold in the laboratory frame of U n , the ionisa- 
tion potential of the target atom in the nth level. This leads 
to a threshold energy in the centre-of-mass frame, and thus a 
lower limit of integration in the calculation of the rate coef- 
ficient, of E m i n - U„/2. The threshold in the centre-of-mass 
frame should be the ionisation energy U„, and in the labo- 
ratory frame should be 2U„. The laboratory frame threshold 
is larger since half of the kinetic energy is, due to conser- 
vation of momentum, unavailable, since it is locked in the 
centre-of-mass motion. 

Furthermore, we experienced difficulties in reproducing the 
results for the rate coefficients. While we were able to per- 
fectly reproduce the cross sections of Soon from the expres- 
sions, when we calculated rate coefficients by direct numer- 
ical integration over the Maxwellian distribution of veloc- 
ities, we encountered problems. Following exactly the for- 
mulae presented in Soon's paper, for the cross sections and 
the rate coefficient integrations, (i.e. making no attempt to 
correct the apparent problem with the thresholds, thus adopt- 
ing E = E\- d b and E m [ n = U„/2) we found significant differ- 
ences between our results and the directly integrated results 
plotted in Soon for the rate coefficients at low temperature 
(T S 10 5 K) and low principal quantum number (n ^ 10). 
However, the results are in good agreement for high T and 
high n. Similar discrepancies were seen with the results from 
Soon's analytical expressions for the rate coefficients. To in- 
vestigate further, we analytically integrated the low energy 
cross section expression for ionisation (Soon's eqn. 10), this 
regime dominating at the temperatures of interest, with the 
aid of the symbolic computing package Mathematica. The 
results were in perfect agreement with our numerical results. 
At the temperatures of interest in this work (T ~ 5000 K), 
the differences between the results are significant, with our 
rate coefficients often 1-2 orders of magnitude smaller than 
those of Soon in this temperature regime. We note that this 
problem cannot be explained by instead adopting E as the 
kinetic energy in the centre-of-mass frame E cm , or a sim- 
ple correction of the threshold in the integrations, as both 
these changes further reduce the rate coefficients at the tem- 
peratures of interest, since the threshold is shifted to higher 
energies. 

Mihajlov et al. /[T996\\2004\\2005i : For the excitation from lev- 
els 4 < n < 10 where n' — n < 5, rate coefficients from 
Mihajlov et al. (120041 120051 ) are available. For the ionisa- 
tion of levels 4 < n < 10, rate coefRcents are available from 
Mihajlov et al. (1996). These semi-classical calculations are 
based on the quasi-resonant energy transfer mechanism (e.g. 
Janev & Mihajlov 1979) and account for processes (U-|7]i, 
including those involving electron transfer. 

Comparison of the various data sets, where they overlap, is 
instructive. We begin with the case of ionisation from the ground 
level, as this also illustrates the problems with the rate coeffi- 
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cients from Soon. A comparison of the rate coefficients from 
different sources is given in Fig. [2] First, we note the appar- 
ent problems with the Soon cross sections and rate coefficients. 
At T S= 10 5 K, our numerically integrated rate coefficients are 
typically about an order of magnitude smaller than those given 
by Soon's analytic approximation, and those numerically inte- 
grated by Soon (not shown, see Fig. 1 1 of Soon's paper), ft is 
also demonstrated that adopting E = E cm worsens the disagree- 
ment. We note also that at the temperatures of interest for the 
photosphere (10 3 < T < 10 4 K), the Drawin formula gives much 
lower rate coefficients than Soon. 

Comparison of the cross sections with experimental results 
available in the energy range E » 10 2 — 10 s eV (see Fig. 1 of 
Soon [1991 and Fig. 3 of Kunc & Soon [T99TT ), shows the Soon 
cross sections are in much better agreement with experiments 
than those of Drawin. This conclusion would not be greatly 
affected by assuming E = E cm , even though the comparison 
with experiment would be worsened for the Soon cross sections. 
However, the rate coefficients at the temperatures of interest are 
severely affected by such a change, since the cross sections near 
the threshold (E cm = E^/2 S 13.6 eV) dominate. In fact, the 
Soon rate coefficients (E = £i a b) are only larger than those of 
Drawin due to the incorrect placement of the threshold. If we 
correct the threshold of the Soon cross sections by assuming 
E = E cm , then the rate coefficients become even smaller than 
those of Drawin (see Fig. |2|. We note, it is dangerous to infer too 
much about the accuracy of the cross sections near the threshold 
from the comparison with experiments at much higher energies. 

Comparisons for ionisation from some selected excited 
states are plotted in the upper panel of Fig. [3] This has 
the advantage of enabling comparison with the results of 
Mihajlov et al. d 19961 1, which we judge to be much more reli- 
able than the classical theories. An immediately noticeable re- 
sult is that for the temperatures of interest, the various sources 
do not follow the same relative behaviour as for ionisation from 
the ground state. For excited states, it is clear that the Drawin 
recipe is typically several orders of magnitude larger than the 
Soon results, irrespective of chosen integration recipe. In gen- 
eral (except for the low-lying states), the Mihajlov et al. results 
lie somewhere in between, smaller than Drawin, yet larger than 
Soon. 

The results for excitation, plotted in the lower panel of Fig.[3j 
are similar. One may also compare the rate coefficients for 
n — 2 — > 3 computed from the Bates & Lewis ( [19551 1 data via 
the detailed balance relation, with those from the general classi- 
cal recipes. However, one must realise that the comparison is of 
limited value since the curve crossing mechanism at work here 
is not reflected in the classical models. Any agreement would be 
purely coincidental. If we compare results at 8000 K, the Drawin 
formula gives a rate coefRcent over an order of magnitude larger 
than the Bates & Lewis value. The analytical formula of Soon 
gives a value a factor of two smaller, while numerical integra- 
tions with E = E\. d \, and E = E cm are roughly one and three 
orders of magnitude smaller, respectively. 

These comparisons clearly demonstrate the present uncer- 
tainties on the rate coefficients due to hydrogen collisions in 
many cases. It is our judgement that the rate coefficients of 
Mihajlov et al. ( [19961 [20041 120051 ), and Bates & Lewis ( [19551 ), 
are to be preferred over the two classical recipes. Thus, these 
data are adopted wherever possible, namely for the ionisa- 
tion from 4 < n < 10 and excitation for the same levels where 
n' — n < 5, and for the case of the n — 2 — > 3 transition. 
Unfortunately, given the discussed problems with the two classi- 
cal recipes, we have no real reason to prefer one over the others. 
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Fig. 2. Comparison of rate coefficients for collisional ionisa- 
tion of H from the ground level due to ground state H atoms, 
H( 1 s) + H( 1 s) -» H( 1 s) + H + + e. The full line is our resu lt from 
direct numerical integration of the cross sections of Soon ( 1992 ), 
where E = E\^. The dot-dashed line shows the result if we in- 
stead assume E is the kinetic energy in the centre-of-mass frame 
E = E cm , rather than in the laboratory frame. The squares show 
results from analytic integration of Soon's eqn. 10, the low en- 
ergy cross section expression. The dashed line shows the result 
from Soon's analytic expression for the ionisation rate coeffi- 
cient. The dotted line shows the result from the analytic expres- 
sion of Drawin ([T968l[l969]). 

For this reason we have chosen to perform calculations for four 
models. We define the following four recipes for calculating the 
remainder of the required rate coefficients: 

HD The Drawin rate coefficients are used. 
HSa Soon's analytical expressions (his equations 15 and 16) are 
used. 

HSn Rate coefficients obtained from numerical integration of 
Soon's cross section expressions, assuming E - E^, are 
used. 

HSnc Rate coefficients obtained from numerical integration of 
Soon's cross section expressions, assuming E = E cm , are 
used. 

We will refer to these recipes in what follows. Our main aim will 
be to determine how the uncertainties in the atomic data translate 
to uncertainties in the non-LTE modelling. While we hope this 
may give us some guidance, we should point out that there is no 
reason to believe that these recipes bracket the true values. 

Mutual neutralisation 

H + +H ^±H(n = 2,3) + H(ls) (8) 
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H(n) + H(1s)^H* + e(n=4,7,10) 



3000 4000 5000 6000 7000 8000 9000 10000 
T[K] 




3000 4000 5000 6000 7000 8000 9000 10000 



Table 1. Polynomial fits to the rate coefficient for collision pro- 
cesses. The rate coefficients (crv) are given by equation[TJin units 
of cm 3 s _1 for the listed process. The fits were made to results 
over the range 1000 < T < 20000 K. 



n or nl 


a\ 02 «3 


(34 




H(nZ) + e -> H + + 2e 




Is 


-2324.8 +1642.6 -391.96 


+31.41 


2s 


-582.7 +409.9 -98.28 


+7.92 




H(n = 2) + H(ls) -> H(n = 3) + H(ls) 


2 


-341.5 +233.9 -55.46 


+4.388 




H + +H" -> H(n) + H(ls) 




2 


-19.08 +8.450 -2.500 


+0.249 


3 


-6.074 +0.048 -0.222 


+0.031 



2.1.3. LTE model atom 

It will be of interest to compare the results of our non-LTE calcu- 
lations with the LTE results. To consistently simulate LTE con- 
ditions, LTE results are calculated using the same codes using 
a model atom with extremely large electron collision rates such 
that all departure coefficients are unity. 



Fig. 3. Comparison of collisional ionisation (upper panel) 
and excitation (lower panel) rate coefficients due to ground 
state H atoms for selected excited states. The data from 
Mihajlov et al. (I1996I l2004i l are plotted as thick full lines. The 
remaining lines follow the scheme used in Fig. [2] i.e., the nu- 
merically integrated data from Soon (1992, E — £) a b) are plot- 
ted as normal full lines. The dot-dashed line shows the result 
if we instead assume E = E cm . The data from Soon's ana- 
lytic expressions are plotted as dashed lines, and the data from 
Drawin ( 1968 1969) are plotted as dotted lines. The dot-dashed 
line for ionisation from n — 7 is under the full line for n — 10. 

This charge exchange process was included employing the cross 
sections of Fussen & Kubach (1986), which have been com- 
puted using a quantum close coupling treatment. This process 
may occur only into the n = 2,3,4 levels as only these states 
have avoided crossings with the ionic state potential. As dis- 
cussed above, the n — 4 crossing is at large internuclear dis- 
tance, and is thus inefficient. The n — 3 level, in particular, has 
a significant charge exchange cross section. The rate coefficients 
for the mutual neutralisation were determined from these cross 
sections, and fits to these results are given in Table Q] Note, as 
H~ is not included in our model atom, in calculating the colli- 
sion rates for this process we have assumed LTE populations for 
FT. As this process could only conceivably be important in the 
deep photosphere where PL is in LTE (see discussion in ® this 
approximation should be adequate. 

Penning ionisation involving H(n=2) 

H(« = 2) + H(n = 2)^PI(ls) + H + +e. (9) 

The quantal calculations of Bates et al. d 19671 1 were employed, 
specifically (crv) = 2.1 x 10~ 8 r~^ 6 cm 3 s _1 for the forward 
reaction. Such collisions are expected to be rare, but we have 
modelled this process anyway as it is easily included. 



2.2. Model stellar atmospheres 

Statistical equilibrium calculations as described above were car- 
ried out for a semi-empirical solar atmosphere model including 
a chromosphere, and a number of theoretical stellar atmosphere 
models. Calculations were performed for the semi-empirical 
model of the quiet Sun of Maltby et al. (1986), hereafter re- 
ferred to as the MACKKL model. The theoretical models are 
computed under the assumptions of ID plane-parallel geometry 
and LTE using the MARCS code (Gustafsson et al. 119751 Asplund 
et al. I1997I I. Scaled solar abundances are used based on solar 
values from Grevesse & Sauval (1998), except for C which is 
taken from Allende Prieto et al. (2002). Convection is treated 
under the mixing-length theory with the parameters a = 1.5 
and y = 3/4jt 2 , which are the default values for MARCS. See, 
for example, Henyey et al. {1965 ) for details. As mentioned ear- 
lier, various authors have suggested revisions of these parame- 
ters, particularly a, to better reproduce Balmer line profiles in 
solar-type stars. However, this investigation is a study of line 
formation processes and this is not important here. 

3. Results 

We begin our study of Balmer line formation with the case of 
the Sun using the MACKKL semi-empirical solar model with 
chromosphere. This permits us to investigate the effect of vari- 
ous collision processes in both the photosphere and the chromo- 
sphere, where the wings and line cores, respectively, predomi- 
nantly form. It is worth emphasising that we study line formation 
within a given model atmosphere, and any feedback of our dif- 
ferent model atom on the semi-empirical temperature structure is 
neglected (see, e.g., Fontenla et al. 2006). Finally, we investigate 
how the results for the line wings translate to other late-type stars 
using theoretical model atmospheres with no chromospheres. 

In order to gauge the effect and importance of the various 
collision processes, we performed calculations separating the 
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various collision processes. As collisions with electrons are well 
understood to be important, we performed a calculation includ- 
ing only electron collisions. Subsequent calculations were then 
performed including additional collision processes. First, mutual 
neutralisation and Penning ionisation processes were added. As 
we will see, these processes have negligible effects on the results. 
We then added collisions with hydrogen atoms, following one of 
the four recipes discussed above, labelled HD, HSa, HSn, and 
HSnc. For brevity, in the following we will refer to these models 
by these labels, noting that it is always implied that electron col- 
lisions, mutual neutralisation, and Penning ionisation processes 
are included. 

3.1. The solar atmosphere 

Figure|4]presents the departure coefficient results from these cal- 
culations for the MACKKL model. To aid in interpretation of 
these results, Fig. [5] plots the total collision rates for selected 
important transitions. First, we point out the result mentioned 
above, the models including electrons only, and electrons plus 
mutual neutralisation and Penning ionisation processes, are es- 
sentially identical. Penning ionisation involving H(n = 2), pro- 
cess (O, was found to have practically no effect on the statistical 
equilibrium, as would be expected due to the rarity of such col- 
lisions. The mutual neutralisation charge exchange process ®, 
couples the n — 2 and n — 3 states slightly more strongly to the 
ionised state, but has only barely noticeable effects. Secondly, 
we notice that the results using the HD and HSnc models are 
not significantly different from those with electrons only. Finally, 
we see that the results from the HSa and HSn models do show 
significant differences from those where only electrons are in- 
cluded. As such, in the following discussions we focus predom- 
inant on the HSa and HSn models since they show the most 
significant changes; however, this should not be interpreted as 
favouring these models. 

In the photosphere, below the temperature minimum (log t <; 
-3), the « = 1,2 states attain LTE populations even with only 
electron collisions included, but the higher states and the ionised 
state show overpopulation in the upper photosphere. To show the 
situation more clearly, Fig.|6]shows the departure coefficients for 
the photosphere for the n — 3 state, which is representative of the 
effect in the higher states. The departures practically disappear 
when excitation and ionisation by collisions with ground state 
hydrogen atoms, processes (|4]-0, are included following recipes 
HSn or HSa. In fact, test calculations showed that the overpopu- 
lations of the H + and n > 3 states in this region of the atmosphere 
are removed by inclusion of only the ionisation by H atom col- 
lision of the n — 1 state. As can be seen in Fig. [5] the collision 
rate for this process increases by several orders of magnitude, 
and thus the n — 1 and H + states are brought close to detailed 
balance bringing H + to LTE, with the effect trickling down to the 
upper states via their already strong coupling to the continuum. 

We now turn our attention to the effect on the Balmer line 
profiles. It is useful for our discussion to have a quantitative idea 
of where the lines form, and contribution functions for the core, 
core- wing transition and wing for the lower three Balmer lines 
have been plotted in Fig.|7]for reference. 

First, we focus our attention on the line wings which form 
deep in the photosphere at log tsoo > -2 (see Fig. [7]). In Fig. [8] 
the differences between the LTE profiles and non-LTE profiles 
in the line wings for Ha, H/3, and Hy are shown. First, we have 
the notable result, that the lines are not well approximated by 
LTE if only electron collisions are included. Secondly, we see 
that even with hydrogen collisions added, the result depends on 



8000 










7000 










H 6000 










5000 










4000 










10 


n=1 - 


1 




10 


\- 






_ 

n = 2 " 


1 


















n = 3 " 


1.0 






e = e + PI + MN 

. HSq 




0. 1 






. HSnc 

. . . HSn 
HD 




' .0 








n = 4 " 










0.1 










1.0 








H* " 










0.1 











-6 -4 -2 

,0 9 T 500 



Fig. 4. Results for the departure coefficients b of the four lowest 
levels of H and for H + with the logarithm of continuum opti- 
cal depth at 500 nm, log T500, from calculations including differ- 
ent collision processes for the MACKKL semi-empirical solar 
model. In the upper panel the temperature stratification is shown 
for reference. The full line shows the result where only electron 
collisions are included. The difference if Penning ionisation (PI) 
and mutual neutralisation (MN) processes are included, in addi- 
tion to the electron collisions, would not be noticeable on this 
plot. The additional lines show the results if all collision pro- 
cesses are included, where the various different recipes for the 
hydrogen collisions are used. 

which recipe for the rate coefficients is used. The results for the 
HD and HSnc models are very similar, and show only a small 
difference from the case with only electrons. For the HSn model, 
the line profiles are close to the LTE profiles in the wings, but 
still some small differences exist. For the HSa model, the line 
wings are in good agreement with those predicted in LTE, and 
are indeed formed in LTE as the departure coefficients and S v /B v 
are unity in the region where the line wings form. Higher Balmer 
lines show similar results. 

In Fig. [8] it is notable that even in the case of the HSa model, 
which attains LTE to high accuracy in the line wings, in the inner 
wing region of Ha around 1 S A/i S 2 A there are very small 
(^ 0.1 %) departures from the LTE profile. This is due to small 
contributions from the chromosphere to this part of the line, as 
demonstrated in Fig. [7] showing the small contributions from 
AA = 0.5. This is shown more clearly in Fig. [9] which compares 
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Fig. 5. Results for the total collision rates between selected lev- 
els with the logarithm of continuum optical depth at 500 nm, 
log rgoo, from calculations including different collision processes 
for the MACKKL semi-empirical solar model. The full line 
shows the result where only electrons are included. The dashed 
and dash-dot lines show the results from the HSa and HSn mod- 
els, respectively. The HSnc and HD models are not shown as 
they show only minor differences from the model including only 
electron collisions. Penning ionisation and mutual neutralisation 
processes do not affect these transitions. Note, the n — 2 — » 3 
transition is in fact calculated using the Bates & Lewis data in 
all cases. The small differences arise due the differences in pop- 
ulations. 
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Fig. 6. Closer view of the departure coefficients b of the n — 3 
state with the logarithm of continuum optical depth at 500 nm, 
log r 5 oo, from calculations including different collision processes 
for the MACKKL semi-empirical solar model, in the photo- 
sphere. Lines have the same meanings as in Fig. [4] 



Fig. 7. Contribution functions to the flux with the logarithm of 
continuum optical depth at 500 nm, log T500, for Ha, H/3, and 
Hy, at three detunings AA ~ 0, 0.5, 5 A chosen to show the be- 
haviours in the core, core-wing transition region, and wing re- 
spectively. These were calculated using the HD model, but do 
not significantly depend on this choice. 



the calculation for Ha using the full MACKKL model with chro- 
mosphere, with a calculation employing the MACKKL model 
where layers above the temperature minimum (logTsoo < -3.6) 
are removed. The far wings show excellent agreement while 
small differences are see in the inner wings. The question if 
there is an effect of the chromosphere on the formation of the 
line wings is of importance. The use of the Balmer line wings as 
effective temperature indicators depends on the fact that they are 
formed in the photosphere, which can be modelled theoretically 
from the basic stellar parameters (usually T e g, log g and chemi- 
cal composition). While the effect seen here in the solar case is 
extremely small, it does perhaps warn of possible effects on the 
inner wing region in more active stars. As an aside, we note that 
Przybilla & Butler (2004b) found that the Ha line profiles, in- 
cluding the cores, were almost identical when computed in this 
manner (MACKKL with and without chromosphere), though 
with quite different line formation physics in each case. We re- 
produced this result if only electron collisions are included (or 
models HD or HSnc); however, as seen in Fig. [9] this result does 
not occur if the HSa model is used (or HSn). 

We now turn our attention to the line core. As shown in 
Fig. |2l the lower Balmer line cores form predominantly in the 
chromosphere, where deviations from LTE are expected as the 
density is low and collisions are inefficient in thermalising the 
gas. We note that the assumption of a static model and statistical 
equilibrium are known to be incorrect in the chromosphere (e.g. 
Carlsson & Stein l2002t Leenaarts & Wedermeyer 2006, and ref- 
erences therein). However, such a calculation may indicate if hy- 
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Fig. 8. The effect of different collision processes on the wings of 
the lower Balmer lines in the semi-empirical MACKKL model 
for the solar atmosphere. The difference between the non-LTE 
line profile Fnlte and the LTE profile Flte is plotted in units of 
per cent of continuum flux. The full lines show the results when 
only electron collisions are included. The remaining lines show 
the results when all collisional processes are included, where 
various different recipes for the hydrogen collisions are used. 




Fig. 9. Comparison of Ha calculated for the MACKKL model 
both with and without chromosphere. The calculations were 
done for the HSa model. 



drogen collision processes could conceivably be important in the 
chromosphere, and in the formation of the line core. Compared 
with when electrons only are considered, the HD and HSnc mod- 
els have basically no effect on the line core. When the HSa or 
HSn models are used the effect on the predicted line cores is 
small, of the order of a few per cent for Ha and H/3. Figure [10] 
shows that the effect on the predicted line cores for Ha and H/3 
for the case of the HSa model; the HSn results are similar. 



Fig. 10. Results for the Ha and H/3 line cores in the MACKKL 
model, with electrons only or the HSa recipe. 

3.2. Late-type stellar atmospheres 

We now investigate how these results for the solar atmosphere, 
particularly for the Balmer line wings, translate to other late-type 
stars, especially the cases of metal-poor and giant stars where 
collisions may be less efficient. Calculations were performed for 
six theoretical models covering a range of different late-type stel- 
lar atmospheres, which are summarised in Table [2] The results 
for the Ha profile are given in Fig. QT| For reference, we in- 
cluded a MARCS solar model, and the results for the line wings 
are seen to be similar to those for the MACKKL model. We 
also calculated for solar metallicity giant, turnoff star, and super- 
giant models. The results are of a similar form to the solar case, 
although the magnitude of the differences from the LTE case 
for the electrons only, HD, HSn, and HSnc models increases. 
The HSa model attains LTE in the line wings in all three cases. 
Calculations were also performed for two metal-poor models. 
The first is an extremely metal-poor giant star corresponding to 
the case of HE 0107-5240, and the second a metal-poor turnoff 
star corresponding to the typical example HD 84937. The results 
are again similar, with the magnitude of the differences from the 
LTE case for the models, where they exist, being different to 
the corresponding solar metallicity case. The electrons only and 
HD model cases show increased differences from LTE, while the 
HSn model results show decreased differences. Once again, the 
HSa models attain LTE in the line wings. 

4. Discussion 

As mentioned in Sect. [1] in most past work on the statistical 
equilibrium of hydrogen in stellar atmospheres, only collisional 
processes involving electrons have been included. This is based 
on the fact that at the thermal velocities of interest electron colli- 
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Table 2. Details for theoretical models used. Unless otherwise stated, models used scaled solar abundances as described in the text. 



Model 






[Fe/H] 


notes 




[K] 


[dex] 


[dex] 




Solar model 


5777 


4.44 


0.0 




Solar metallicity giant 


5100 


2.20 


0.0 




Solar metallicity tumoff star 


6340 


4.00 


0.0 




Solar metallicity F supergiant 


6000 


1.50 


0.0 




Very metal-poor giant 


5100 


2.20 


-5.3 


a-enhanced (+0.4 dex), [C/Fe] = 4.0, [N/Fe] = [O/Fe] = 2.3 


Metal-poor turnoff star 


6340 


4.00 


-2.3 


cf-enhanced (+0.4 dex) 



sions are expected to be non-adiabatic and have substantial cross 
sections, while collisions with heavier particles are expected to 
be nearly adiabatic and thus have small cross sections. In the 
photosphere of solar-type stars, neutral hydrogen atoms typi- 
cally outnumber electrons by a factor of 10 4 , and can be orders 
of magnitude greater still in metal-poor stars. For example, for 
the HSa and HSn models, the rate coefficients at 5000 K for the 
n = 1 — > 2 transition due to electrons are only of the order 15 
and 2500 times greater than those due to hydrogen collisions, re- 
spectively, and thus hydrogen collisions may even be dominant 
in certain regions as seen in Fig. [5] For the rate coefficient for the 
n — 2 — > 3 transition the difference is of the order 10 6 , and thus 
hydrogen collisions start to become important for this transition 
in metal-poor stars. Thus, sheer weight of numbers could lead 
to the importance of hydrogen atoms collisions, which has been 
seen in the presented calculations. 

It is appropriate to note here, that the above arguments speak 
against the possibility of inelastic collisions with protons being 
of any significance. As collisions with protons are expected to be 
adiabatic and the proton abundance will be less than the electron 
abundance, such processes should always be insignificant with 
respect to electron collisions. However, one must also note the 
limited validity of the Massey criterion. The adiabatic approx- 
imation breaks down in the case of small energy separations, 
such as occurs at avoided crossings of potential energy curves 
(e.g. Bohm [T95~n page 505). 

We have shown that mutual neutralisation processes and 
Penning ionisation are unimportant in the statistical equilib- 
rium of hydrogen in late-type atmospheres. Unfortunately, from 
the results presented here the role of collisional excitation 
and ionisation by hydrogen is unclear. Two models, HD and 
HSnc, predict that hydrogen collisions are relatively unimpor- 
tant. However, the HSn and HSa models predict that hydrogen 
collisions are important, and in the latter case predict that Balmer 
line wings form in LTE conditions. It should be noted that the lat- 
ter two models have severe problems in their formulation; in par- 
ticular, they are based on cross sections which have physically 
incorrect thresholds. On the other hand, the HD model does not 
compare well with experiment at high energies, and the HSnc is 
essentially a correction to the HSn model with no other physical 
justification than to obtain the correct threshold. It is therefore 
urgent that more physically sound theoretical results or exper- 
imental data are obtained for these processes at near threshold 
energies. Excitation and, particularly, ionisation from the lowest 
lying states are most vital. 

Naturally, these results cast doubt on the usefulness of the 
Balmer line wings as reliable temperature indicators for late- 
type stars. The presumed advantage of Balmer lines has been 
that they are spectroscopic, and therefore unaffected by redden- 
ing (as opposed to photometry), and formed in LTE (as opposed 
to, e.g., Fe I excitation equilibrium). If the assumption of LTE is 
invalid, then without high quality atomic collision data, so that 



the departures from LTE can be accurately predicted, the tem- 
peratures derived from Balmer lines will be subject to errors. 
The region of the Balmer line most sensitive to effective tem- 
perature, and that usually used and most important in fitting the 
lines, is the region around 2 S A/l S 6 A. As a rough rule of 
thumb, a one per cent change in the inner wings of Ha corre- 
sponds to 60-100 K, the exact number depending on stellar pa- 
rameters (Fuhrmann et al. |1993t Barklem et al. 2002). For the 
cases considered here, we see that the possible departures from 
LTE (given by models HSnc or HD) in this region are of the 
order of a couple of per cent, thus leading to changes (raising) 
with respect to the temperatures determined in the LTE approxi- 
mation of order 100 K. 

It is worth considering how such changes would impact com- 
parison with observations. Such corrections can probably not 
be confirmed or ruled out on the basis of observational con- 
straints, in particular, agreement with other temperature indica- 
tors (e.g., Barklem et al. 2002), given the errors in the analyses 
are of similar magnitude. A slight weakening of the predicted 
solar Balmer line profiles would improve agreement with obser- 
vations in terms of line strength, where a temperature of order 
50 K too low is found using the LTE approximation for the case 
of the MARCS models and our adopted broadening theory (see 
Barklem et al. 2002). However, Fuhrmann et al. ( 119931 ), using 
different model atmospheres and broadening theory, have found 
temperatures in agreement with the known solar value. This em- 
phasises the difficulty in testing the atomic data via comparison 
with observations, due to the uncertainties in the model atmo- 
sphere and line broadening physics. Even if semi-empirical at- 
mospheres are used, the uncertainties due to assumptions such as 
ID stratification and uncertainties in the line broadening physics 
remain. Furthermore, it is somewhat inconsistent to compute for 
a semiempirical model that has been calculated using a different 
model atom, and this could possibly result in minor errors (see 
Fontenla et al. 2006). Regarding line shape, since the very in- 
nermost parts of the wings are most significantly affected, the 
shapes of the lines would also be changed. This would cer- 
tainly lead to discrepancies with observed line shapes, though 
this could probably be resolved with a revision of mixing-length 
convection parameters, at least for theoretical models. Any ob- 
servational tests of the atomic data must be treated with care, due 
to the assumptions and additional uncertainties such those of a 
ID static atmosphere (often modelled in LTE) and line broad- 
ening physics. While detailed and careful comparisons of the 
different atomic model predictions with observations, preferably 
using 3D models, would indeed be interesting, we believe that 
the preferred approach to dealing with this problem is to obtain 
improved theoretical and/or experimental data for the hydrogen 
collision processes, as discussed above. 
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Fig. 11. The effect of different collision models on the wings of 
Ha for various MARCS models. The difference between the non- 
LTE line profile Fnlte and the LTE profile Flte is plotted in 
units of per cent of continuum flux. The lines have the same 
meanings as in Fig. [8] but the long-dashed line (HSNc model) 
is not plotted, as it is practically coincident with the dotted line 
(HD model) in all cases. Each plot is labelled by the model pa- 
rameters r ei f/logg/[Fe/H]. 
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